PROCEEDINGS 

°F SCIENCE 



mb and /b, from a combination of HQET and QCD 



o 
o 

(N 

CD 

00 
(N 



> 

in 
\o 
o 

o 

o 

CD 



X 



Damiano Guazzini*, Rainer Sommer 

DESY, Platanenallee 6, D- 1 5738 Zeuthen, Germany 



E-mail: Damiano.Guazzini@desy.de, Rainer.Sommer@desy.de 

Nazario Tantalo 

INFN sezione "Tor Vergata", Via della Ricerca Scientifica 1, 1-00133 Rome, Italy & 
Centro E. Fermi, Compendio Viminale, 1-00184 Rome, Italy 



E-mail: Nazario . Tantalo@roma2 . inf n . it 



We compute the mass of the b-quark and the Bg meson decay constant in quenched lattice QCD 
using a combination of HQET and the standard relativistic QCD Lagrangian. We start from a 
small volume, where one can directly deal with the b-quark, and compute the evolution to a big 
volume, where the finite size effects are negligible through step scaling functions which give the 
change of the observables when L is changed to 2L. In all steps we extrapolate to the continuum 
limit, separately in HQET and in QCD for masses below m^,. The point m\, is then reached by an 
interpolation of the continuum results. With ro = 0.5 fm and the experimental and K masses 
we find /b^ = 191(6) MeV and the renormalization group invariant mass Mb = 6.89(11) GeV, 
translating into mb(mb) = 4.42(7) GeV in the MS scheme. 
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1. Introduction 



The b-quark mass is a fundamental parameter of QCD, and its accurate knowledge is needed 
for theoretical predictions of B meson decay rates. The understanding of the latter is a very active 
field of high energy physics research. At the same time the B meson decay constant plays a crucial 
role in the description of these phenomena. 

We focus our attention on the pseudoscalar Bg meson, a system characterized by two different 
scales: the heavy quark mass (mb ~ 5 GeV) and the typical QCD scale. The mass of the strange 
quark is around or below the latter. We fix it to its physical value through wjk as in 



2. HQET and step scaling method (SSM) 

We deal with these two scales in (quenched) lattice QCD with the SSM introduced in [§, |3|], 
but constraining the large mass behaviour by HQET [Q]. The computation of an observable 0{m\^) 
using the SSM is based on the identity 

OK,L.) = OK,Lo) • • • ' (2.1) 

where wih stands generically for a heavy quark mass whose precise definition is needed only later. In 
order to be able to extract each factor in the continuum limit, the starting volume Lq has to be small 
enough to properly account for the dynamics of the b-quark, using a relativistic 0(a) -improved 
action. A good choice is Lq = 0.4 fm [^, ^, where easily lattice spacings of a » 0.012 fm can be 
used. (Physical units are set using = 0.5 fm |6[ ^). Furthermore, Leo has to be large enough 
such that finite size effects in 0{mh,Lco) are negligible. In practise we will use Loo ~ L2 = 1.6 fm. 
We will choose a fixed ratio s = L,/L,_i in the step scaling functions 

ao(mh,L,-) = ^(^^^^._^). (2.2) 

The number N and the scale ratio s of the steps are in principle dependent on the considered 
observable and on the desired level of accuracy. It has been seen [^, |3|] that (A'^,^) = (2,2) is a 
suitable choice for the mass and decay constant of the Bs meson. 
In HQET the step scaling functions are expanded as 

ao(mh,L,) = a(°)(L,) + + (-1 ) (2.3) 

at fixed L,. We wiU see that the correction terms to the leading order are small for the masses of 
interest. 

We first consider the case of a finite volume pseudoscalar meson mass, C?(mh,L) = Mps(/Mh,L), 
which will be defined in the following section. In this case, a^°' = 1 and the first non-trivial term 
O^Q^ is computable in the static approximation of HQET. We further define 
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as the natural non-perturbative dimensionless mass variable. The step scaling function for the 
meson mass is then written as 

G^{x,U) ^ =\ + GT{Li)-x + 0{x^), x = x{m^,Lt). (2.5) 

It is defined for all x,L. The idea for its numerical evaluation is to compute a^^{L) explicitly in the 
static approximation and fix the small remainder by the relativistic QCD data with quarks of masses 
of the physical charm quark and higher. In other words we interpolate to the physical b-quark mas. 
With the experimental mass of the Bj meson, Mq^ = 5.3675(18) GeV we fix X2 = 1 /L2Mq^ and the 
physical points corresponding to the b-quark are then given by 

X2 = 1/(L2MbJ , x,_i = 2am(x,-,L,) -Xi . (2.6) 

The numerical results will have to be evaluated at these points. In the smallest volume we relate 
the meson mass to the renormalization group invariant (RGI) quark mass, Mh, defining 

p(x,Lo) = ^5^Psi^ = p(o)(z.,,)+p(i)(Lo).. + 0(x^). (2.7) 
We thus have the connection of the Bg meson mass and the RGI b-quark mass 

Mb = I J N 7^^—\ 1 TT • (2-8) 

p(x,Lo) • am(xi,Li) • am(X2,L2) 

For the decay constant the step scaling function 



af(.,L,.) ^ , /psK,L,.)^^s feL4^ ^ or{U) + oi'\u) .x + 0{x') (2.9) 

yields straightforwardly the connection between the finite volume decay constant and the infinite 
volume one. Note that the only approximation made in the above equations is to neglect finite size 
effects on mass and decay constant in the volume of linear extent L2. 

3. Finite volume observables 
3.1 Relativistic QCD 

Suitable finite volume observables are defined in the QCD Schrodinger functional [^, ^ with 
a space-time topology xT, where T = 2L and C = C = is chosen for the boundary gauge 
fields, and = for the phase in the spatial quark boundary conditions. 

The 0(a)-improved correlation functions /A(mh,L,xo),/p(mh,L,xo) and f\{mh,L) are defined and 
renormaUzed as in allowing to compute the pseudoscalar meson decay constant 

/psK,L) = /aK'^'^) -^="Vb.(L)"=°°/b. (3.1) 



and the pseudoscalar meson mass 



Mps(mh,L) = ^In 



fAimh,L,L-a) 
fA{mh,L,L + a) 



■Mj,SL)=Mb, (3.2) 
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For all observables computed in relativistic (quenched) QCD we employ the non-perturbatively 
0(a)-improved Wilson action[l^, 11]. The data at finite heavy quark mass were published in [| 



They have been reanalyzed, taking into account the correlation between observables computed on 
the same gauge configurations. The statistical uncertainties on the renormalization constants and 
the lattice spacing are included before performing the continuum limit extrapolations; they do not 
appear as a separate uncertainty. 

3.2 HQET 

In the static approximation of HQET, unrenormalized correlation functions /^'^' and ff^^ are 



defined in complete analogy to the relativistic ones, see [12]. As in this reference, we use the RGI 



static axial current, related to the bare one by a factor Zrgi. It serves to define the RGI ratio , 

FRGi(L)=ZRGi4^ii^, (3.3) 

which is related to the QCD decay constant /ps via 



fvs{mb,L)^DM^s{L) = -2Cps(AMs/Mh) x Frgi(L) + 0(l/mh) . (3.4) 

The function Cps(Aj^/Mh), defined in [p^, can be accurately evaluated in perturbation theory; 
we use the 3-loop anomalous dimension 7^^ computed in [|^]. Just like Zrgi, it is needed only for 
/ps('Mh,Z.o); it cancels out in the step scaling functions. 



In analogy to eq. (|I|) we further define T^i^^{L) = ^\n[ff'\L,L- a) / ff''\L,L + a)] . The 



static step scaling functions then read 

uf'^'iLi) = 1 W^/) ^stat(^.) ^ ^. [p^^^^(^.) _ r,,,,(L;_0] , Li = 2L,-_i . (3.5) 

2 -'RGH^i— 1 ) 



These quantities will be precisely computed by using the static action denoted by HYP2 in [ |14| ] 



(see also [15]), and the corresponding 0(a) -improvement coefficients for the static axial current. 



The regularization independent part of the factor Zrgi is known from []12|], while the regularization 
dependent one is computed in this work. 

4. Numerical results for the b quark mass 

The computation of am{x,L2) is performed at finite quark mass on lattices with j8 = 5.960, 
6.211, 6.420 and resolutions L2/a = 16, 24, 32; the continuum limits for the three heaviest 
quark masses are shown on the left of Figure [l[ For the static step scaling function we took the 



results for L = L2 from an extension! 16] of the work of the ALPHA collaboration[rM], while in the 



intermediate volume (Li) we simulated lattices with 5.960 < j8 < 6.737. The continuum limit 

ar(^2) = 1.549(33), (4.1) 

is used in the interpolation of am{x,L2) between values of x corresponding to about the mass of 
the charm quark and the limit <7m(0,L2) = 1. It constrains the slope of the fitting curve to the cone 
shown in Figure |l|. The result of the quadratic fit in x reads 

am(x2,L2) = 1.0328(11), (4.2) 
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Figure 1: Continuum limit extrapolation and interpolation of (7m('«,i2) 

hardly distinguishable from a purely static result. Analogously the interpolation of the step scaling 
function for the intermediate volume gives 



a„,(^i,Li) = 1.0092(18). 



(4.3) 



In the small volume only the relativistic data are needed to establish a finite volume relation- 
ship between the meson and the heavy quark masses. The renormalization is non-perturbatively 
achieved through the renormalization factor Zyiigo) and the 0(a) -improvement terms computed in 



[118|, [19|, gO|]. Using eq. (|2.8D, the interpolated value 

p(xo,Lo) =0.748(11), 



(4.4) 



is combined with the above step scaling functions to find the scale and scheme independent number 
Mb = 6.888(105) GeV => m^^{m^;^) =AA2\{61)GqW . (4.5) 

5. Numerical results for the decay constant 

For the computation of Of{x,L2) the relativistic data originate from the same gauge configura- 
tions used earlier, while in the static case the decay constant in the bigger volume, 

}rgi(L2) = -4.63(19), (5.1) 

was again computed and extrapolated to the continuum limit as an extension [ |T6| ] of [17]. The 
continuum extrapolation of the same quantity in the intermediate volume (L = L\) is shown on the 
left of Figure ^. The result 

>'RGi(ii) = -1.628(19) (5.2) 
is used together with ( |5.1[ ) and the relativistic data, as shown on Figure ^ (right), to get 

of'-\L2) = 1.006(44), Of{x2,L2) = 0.974(30) . (5.3) 

Similarly, but by extrapolating the step scaling function to the continuum limit rather than Frgi(^i) 
and Frgi(^o) separately, we obtain 



: 0.4337(44), 



af(;t;i,Li) =0.4260(31). 



(5.4) 
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Figure 2: Continuum extrapolation of FsF(i'i ) and interpolation of CTf (x,L2) 



With the small volume results (see Figure ^ 

-/b (Lo)a/l^Mps(Lo) 

Yrgi{Lo) = -1.347(13), Y^^ixoM = 2C^^(j^^j m,,) = -1-280(17) , (5.5) 

we finally arrive at the result 

/b^ = 191(6) MeV. (5.6) 



6. Conclusions 

The combination of the Tor Vergata strategy to compute properties of heavy-light mesons 
[@, with the expansion of all quantities in HQET [Q], changes extrapolations in the former com- 
putations into interpolations. As expected, our numerical results demonstrate that these are very 
well behaved. Indeed the higher order mass dependence of the step scaling functions is very weak, 
and in all but one steps the static approximation alone gives very accurate results. In the one 
exception (yps(xo,Lo), Fig. ^ the 0(1 /wJb) corrections are around 5%. 

Our results do not suffer from any systematic errors apart from the use of the quenched ap- 
proximation; small systematic errors quoted in [^, ||] for the extrapolation uncertainties have been 
eliminated. Our results are in agreement with the ones of [|,|,|l7Hll], within the errors. 



Concerning dynamical fermion computations, the challenge in this strategy is to simulate in a large 
volume (such as L2) with small enough lattice spacings, where quark masses of around mcharm and 
higher can be simulated with confidence. 
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Figure 3: Continuum extrapolation of Ys,f{Lq) and interpolation of the decay constant on the small volume 
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